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Abstract 



We extend and generalize the seminal work of Brandenberger, Huang and 
Zhang on the formation of strings during chiral phase transitions Q and dis- 
cuss the formation of abelian and non-abelian topological strings during such 
transitions in the early Universe and in the high energy heavy-ion collisions. 
Chiral symmetry as well as deconfinement are restored in the core of these 
defects. Formation of a dense network of string defects is likely to play an 
important role in the dynamics following the chiral phase transition. We 
speculate that such a network can give rise to non-azimuthal distribution of 
transverse energy in heavy-ion collisions. 
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I. INTRODUCTION 



QCD with massless quarks possesses U(Nf) R x U(Nj) L symmetry. However, the low 
energy spectrum of QCD does not exhibit this full symmetry and suggests that part of this 
symmetry is spontaneously broken. It is the axial vector part of this symmetry which is 
spontaneously broken at low energy, U(Nf) R x U(Nj)l — > U(Nf)v=R+L = U(Nf)y. Pseu- 
doscalar mesons are the Goldstone bosons of this spontaneous symmetry breaking (SSB). 
The vacuum at low energy/temperature, i.e in the broken phase, is characterized by a non- 
zero value of the chiral condensate. 

In the real world, the U(Nf)A transformations, apart from spontaneous breaking, are 
also explicitly broken because of nonzero quark masses. Due to this the Goldstone bosons 
become massive. Even in the chiral limit, the U(1)a part of U{Nj)a is spontaneously broken 
to Z(Nf) due to anomaly 0. This makes the corresponding Goldstone boson, the rj' meson, 
heavier compared to other pseudoscalar mesons. Finite temperature field theory calculations 
on lattice show that SU(Nf)A is restored above a critical temperature T x in the chiral limit 
and approximately restored for the case of small quark masses. Also lattice results 
show that topological susceptibility sharply drops above T x effectively restoring the U{1)a 
symmetry. Such chiral symmetry restoration could have occurred in the early Universe, 
when the temperature of the universe was above T x . Also, one of the main objectives of 
relativistic heavy-ion collisions is to detect and study such chiral symmetry restoration, or 
spontaneous breaking of the chiral symmetry, in QCD. 

One of the interesting aspects of spontaneous symmetry breaking phase transitions is 
that topological defects may form during these transitions, depending on the structure of 
the vacuum manifold. Topological defects are extended, non-zero energy solutions of the 
equations of motion. Their stability typically originates from a non-zero topological charge. 
There are numerous examples of topological defects in condensed matter systems, such as 
flux tubes in super-conductors, vortices in super-fluids, strings and monopoles in nematic 
liquid crystals etc.. Topological defects like cosmic strings, monopoles and domain walls 
etc. arise in various particle physics models of the early Universe. Even though topological 
defects arise in the symmetry broken phase, in the core of these defects symmetry is restored, 
as if the phase with higher symmetry is trapped in the background of broken phase. 

For the case of small explicit symmetry breaking, these topological field configurations 
are no more static solutions of the theory. Nevertheless they will always form during phase 
transitions, though their interaction dynamics may be dominated by explicit symmetry 
breaking. These defects have structures which are stable against small fluctuations of the 
order parameter (even though the defect may not be a static solution) and can only decay 
by annihilation with defects of opposite topological charge. This makes them long lived and 
can lead to important effects on the dynamics of the system. 

In QCD, in the chiral limit, spontaneous symmetry breaking U(Nf) R x U{Nf)i — > 
Uy(Nf) allows for existence of topological string defects. Formation of topological and 
non-topological string defects during the chiral transition in QCD has been discussed 
by Brandenberger, Huang and Zhang The topological strings in JIJ are associated 
with the SSB of U{1)a of U(Nf)A and the authors argue that they will decay as the 
anomaly effects become substantial. We will show below that the spontaneous breaking 
of U(Nf) R x U(Nf) L — > Uy(Nf) leads to a class of static topological string solutions apart 
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from the topological strings discussed in IJ. Even with anomaly effects which break U(1)a 
of the U(Nj)a further to Z(Nf), the topological strings remain static, although the field 
configuration is different from the case without anomaly effects. In the framework of linear 
sigma model we give the field configurations of these topological string defects in the chiral 
limit including the effects of anomaly Also we will argue that chiral symmetry as well as de- 
confinement are restored in the core of the string which could have interesting implications. 
These topological string solutions become non-static for non-zero quark masses. 

In the next section we briefly discuss the connection between topological defects and 
spontaneous symmetry breaking. In section III we will discuss the essential features of the 
linear sigma model Lagrangian which admits static topological string defects in the chiral 
limit @. Following this in section IV we will give the static string solutions in the chiral 
limit taking parameters from reference [|J. In section V we will discuss how these defects 
will form during the chiral transition in the chiral limit with some remarks for the case of 
non-zero quark masses. In section VI we will discuss the implications of formation of such 
defects. 



II. TOPOLOGICAL DEFECTS AND SSB 

The type of defects formed in a phase transition and the symmetry breaking pattern 
are intimately related. Complete SSB of U(l) symmetry gives rise to vortices and strings 
in 2-D and 3-D physical space respectively. SSB of SU(2) to U{\) gives rise to monopoles 
in 3-D physical space. For the chiral transition one can argue that SSB of chiral symmetry 
U(Nf) x U(Nf) to U(Nf) leads to strings. Energetics of the defects is primarily decided by 
whether the SSB is that of a local or global symmetry. In the case of local symmetry one 
usually has finite energy topological defect solutions, on the other hand for global symmetry, 
like in the chiral transition, energy of the solutions diverges logarithmically with the volume. 

Topological defects arise when the ground state of the theory is degenerate and the 
order parameter space (OPS) consisting of all possible values of the order parameter (OP) 
is topologically non-trivial. For example when £7(1) symmetry is spontaneously broken the 
OP is given by ip = rje 10 with non-zero fixed 77. The OPS in this case consists of all possible 
values of 9 between and 2tc, hence the OPS is a circle S 1 . An OPS has non-trivial n-th 
homotopy group n n if there exist non-trivial mappings from n-sphere S n to the OPS, or 
in other words if one can have non-contractible loops, closed surfaces etc. in the OPS. In 
general topological defects are characterized by the homotopy group it n {OPS) of the OPS. 
Group ir n (OPS) consists of all non-trivial mappings from S n to OPS. Vortices or strings 
are characterized by iri(OPS), monopoles by 7r 2 (0 PS) etc.. For the case of OPS=S' 1 only 
7Ti(S' 1 ) is non-trivial and each element of this group relates to total variation of phase of the 
OP around a vortex or string. The total phase variation around a string or vortex is always 
an integer multiple 2nn of 2n. Because of this, the group 7r 1 (S' 1 ) is isomorphic to the set of 
integers Z. n here is called the winding number which is the topological charge of the string 
or vortex configuration. 

For the SSB of chiral symmetry U(N f ) R x U(N f ) L — ► U(N f ) the OPS topologically is 
equivalent to U(N f ) = SU W* U W . One can show that there exist non-trivial closed loops 



in U(Nf) ||. Closed loops in the OPS can be generated by the Ua,(1) generator. These 
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correspond to topological defects known as the abelian strings . Also there are non-trivial 
loops in OPS generated by linear combination of the Ua(1) and SU(Nf) generators, which 
give rise to non-abelian strings. We consider here the abelian topological strings and will 
discuss non-abelian strings in a separate work. In the next section we will consider the linear 
sigma model Lagrangian |4J] and will give the configuration of the winding one string in the 
chiral limit. 



III. THE MODEL LAGRANGIAN 

The U(Nf)n x U(Nf)i linear sigma model for Nf quark flavors is given by @,[|, 

£($) = Tr (dn&dPQ - m 2 $ f $) - \ x [Tr (V$)] 2 - A 2 Tr (V$) 2 

+ c [Det ($) + Det + Tr [#($ + . (1) 

$ is a complex Nf x Nf matrix field parameterizing the scalar and pseudoscalar mesons, 

$ = T a <p a = T a (a a + m a ) , (2) 

a = 0, ...,Nf — 1 and T a , a ^ 0, are a basis of generators for the SU(Nf) Lie algebra, the 
analogoues of A a for Nf = 3. T = 1 is the generator for the Ua(X) Lie algebra. 
The field $ transforms under the chiral transformation U(Nf)n x U{Nf)i as 

$ — ► U r $ U} , = exp (i T a ) . (3) 

One can rewrite the left-right transformations in terms of vector-axial vector transformations 
with parameters Wy A = (u>" ± u>f)/2. The infinitesimal form of the U(Nf) R x U(Nf) L 
symmetry transformation (^) is 

T a a — > T a a - z u;£ [T a , T b ] 6 + 2 a;^ {T a , T b } 6 . (4) 

$ is a singlet under the U(l)y transformations (exp(iw v T )) . In QCD, it gives rise to a 
conserved charge identified with baryon number. 

In the above model the determinant term takes into account the instanton effect which 
explicitly breaks the U(1)a symmetry 0. It is not clear whether such a term can be 
justified for high temperatures, because vanishing $ makes the instanton effects to disappear 
0, though lattice results show similar behavior for the chiral condensate and topological 
susceptibility . Both drop sharply above the critical temperature T x suggesting that they 
become vanishingly small in a narrow range of temperatures. We expect that for such a 
narrow range of temperatures the basic picture of string defect formation during the chiral 
transition will not be affected. The last term is due to non-zero quark masses, where 
H = T a h a . 

When c = H = A 2 = 0, U(Nf) x U(Nf) is spontaneously broken for m? < and Ai > 0. 
For A 2 7^ U(Nf)A is spontaneously broken to identity. This results in Nf Goldstone 
bosons, for Nf = 3 these Goldstone bosons are the 7r's, i^'s, r] and 7/ . However when just 
c 7^ 0, the U(1)a is further broken to Z(Nf) by the axial anomaly, making the 7/ massive 
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compared to other pseudoscalar mesons. SU(Nf)y is still the symmetry of the Lagrangian. 
All these symmetries are in addition explicitly broken by non-zero quark masses making all 
the Goldstone bosons massive. 

In the absence of anomaly and H = 0, the non-trivial loops in the OPS can be generated 
by exponentiating To or a linear combination of To and T a , a — 1, Nj — 1 0. Formation 
of topological defects, corresponding to the non-trivial loops generated by exponentiating 
To, during the chiral transition have been discussed in For c ^ and H = 0, only the Nf 
equally spaced points %e i2nn/N f , n = 1, 2, N f on the U A (l) circle ($ e id : < 9 < 2tt, $ 
a suitable fixed point in OPS ) are degenerate and all other points have higher effective 
potential or free energy regardless of the values of Aj's. However one can still have static 
string configurations with non-zero topological charge. For H ^ as well, only one point 
on OPS is the absolute minimum of the free energy. There exists no static topological 
configuration in this case, however for small H ^ one can expect topological string field 
configurations to form during the chiral transition. In the following we find the topological 
string solutions for both c = and c 7^ in the chiral limit for Nf = 3. 



IV. THE STRING SOLUTION 

For the abelian topological string solution, it is enough to consider the variation of $ 
generated by To. Topological aspects of the string will not change by considering variation 
of other components of $. However in the core of the string all components of $ do vanish 
in the linear sigma model we are considering here. We assume <3> = T o (0i + ifa), 02 is the 
field for 7/. With this restriction on $, the effective Lagrangian reduces to that of a complex 
scalar field given by the following equation. 

£($) = ^0^0! + i^02^0 2 - V(01, 2 ) (5) 

where 

v(0i,0 2 ) = ^ (0? + 0^ + ^(0? + ^) 2 

- I (0? - 30!02) - 6/i o 0i , (6) 

ho = - l -Tr{HT,). (7) 

We consider here = for which one can solve the field equations for static string 
configurations. The static string configuration is straight and has a certain symmetry around 
the length of the string. For c = the solution of the string has cylindrical symmetry. 
However for c ^ the string does not have cylindrical symmetry, but possesses a symmetry 
of rotation by 27r/3 about the axis of the string. In this case the string is the junction of 3 
domain walls which interpolate between the different Z(3) ground states. In the following 
we consider the string solution first for c = and then for c ^ 0. 

Considering the string along the z-direction, (f>i(x,y) and 02 (x,y) for the static string 
satisfy the following field equations: 
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2 dvfaM 

V 2 2 = gj%M, (8) 



For winding one cylindrically symmetric string solution one considers the following ansatz 
for the string: 



0(f) = 0(r)-. (9) 
r 

Here = J<f>1 + (f>% and r = y/x 2 + y 2 . Now 0(r) satisfies the field equation 

** + i^_4 =m V + A^. (10) 
ar A r ar r z 

We solve the above equation using the parameters ||] corresponding to m a = l.GeV and 
vanishing masses for all pseudoscalar mesons. In the left of FIG.l, we plot — and in the 
right plot the vector field of 0(f). Magnitude is proportional to the length of the vectors 
and the angle these vectors make with the positive x-axis is the U(1)a phase. Configuration 
of U(1)a phase is symmetric around the string in this case. In the core of the string 0(f) 
vanishes restoring the chiral symmetry. 




FIG. 1. Configuration of the string for c = and ho = 0. Plotted is —<ft of the string in the left 
figure and figure on the right gives the vector field plot of the full field <p\{r) + ifoif). 

For the case of c ^ it is difficult to find the static solution by solving the field equa- 
tions. We have found an approximate solution by evolving a cylindrically symmetric string 
configuration by putting an additional dissipative term in the equations derived from (6), 
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FIG. 2. Configuration of the string for c = 2.729 and /io = 0. Figure on the left gives — y + 02 
and figure on the right gives the vector field plot of the full field </>i(r) + ifoir). 

+ V 0i = — 



rf 2 02 . d(j)2 2 9V(0i,0 2 ) 

^ + -V0 2 = ^— , (11) 

where a is the dissipation coefficient. The dissipation term has been used to converge the 
initial configuration (the circular configuration, which is not the correct configuration in 
the presence of the anomaly) to the approximate configuration of string connected with 
domain walls. Once the evolution is done for sufficiently long time the configuration now 
does not evolve with time even after the dissipation is switched off, which suggests that 
the configuration is more or less static and so an approximate solution. On evolving the 
initial cylindrically symmetric solution with the above equations, domain walls connected 
to the string develop. In this case we considered the rf mass to be the experimental value 
m v r ~ 950MeV and set other pseudoscalar meson masses to zero. In the FIG. 2 we show 
the configuration of the string with the domain walls joined to it. In the left figure we 



plot —\J4>i + 4>2, and in the right plot the vector field of <3>(r) = {<j>x{r),<j>^{r))- Again the 
magnitude of the field is proportional to the length of the vectors and the angle these vectors 
make with +ve x-axis is the U(1)a phase. Clearly there is no symmetric distribution of the 
U(1)a phase around the string. Here as well the magnitude of the field vanishes in the core 
of the string restoring the chiral symmetry. 
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V. FORMATION OF STRING NETWORK DURING THE CHIRAL 

TRANSITION 



In the above we have discussed the configuration of static topological string configurations 
in the chiral limit in the linear sigma model. These topological configurations become non- 
static for non-zero quark masses. In this case the phase variation around the string will 
not even have the 3-fold symmetry as in Fig. 2. There will be a domain wall, with phase 
variation of 2im across, connecting a string with topological charge n. Both abelian |T[ 
and non-abelian string defects are expected to form during the chiral transition. In the 
chiral limit one can use the picture of conventional mechanism of defect production, known 
as the Kibble mechanism || to make an estimate of the density of strings. (We mention 
that in situations where the magnitude of the order parameter undergoes huge fluctuations, 
a new mechanism of defect formation becomes effective, as discussed in ||. We do not 
consider such possibilities.) For simplicity we consider here the string configurations which 
are specified by the {7a (1) phase only. However consideration of non-abelian phases will 
increase the probability of formation of defects. 

The basic picture of this mechanism is the following. For a second order phase transition, 
after the transition, the physical space develops a domain like structure, with the typical 
size of the domain being of the order of a relevant correlation length £ (which depends on 
the dynamics of the transition). Inside a given domain, the {7(1) phase is roughly uniform, 
but varies randomly from one domain to the other. The phase in between any two adjacent 
domains is assumed to vary such that the variation of the phase is minimum. This leads 
to the formation of string defects in the junction of three or more domains. Numerical 
simulations [10] show that on an average the expected number of strings per domain is 0.88. 
However for the case of explicit symmetry breaking other mechanisms can also contribute 
to the formation of these defects. In generic situations formation of topological defects 



is enhanced in systems with explicit symmetry breaking pl [. However, this enhancement 



crucially depends on the nature of the dynamics of phase transition. In contrast, Kibble 
mechanism allows for estimates of defect production (per domain) which are reasonably 
independent of the dynamical details. Therefore, for a rough estimate for the density of 
defects one can still use the Kibble mechanism. Thus, we will assume a domain-like structure 
after the phase transition, and assume a random distribution of the U(1)a phase between 
[0, 27r] in the domains. 

So immediately after the chiral transition there will be a network of strings in the broken 
phase. For an infinite system one would expect both open and closed string to form, for 
example during the chiral phase transition in the early universe. But for transition taking 
place in a small volume with outer region in the broken phase, which is the case of heavy-ion 
collisions, one will expect only string loops to form. Interestingly the number of string loops 
close to the boundary will be larger as discussed in ref. [I2"|] . The usual mode of decay of 
these string loops is via shrinking and breaking of bigger loops to smaller loops. 

One of the implications of formation of strings during the chiral transition can be due 
to the interplay between the chiral transition and deconfinement. Recent lattice results 
suggest that chiral condensate and the Polyakov loop are strongly correlated ||14J| . In [15| 
it was proposed that the chiral condensate acts as an effective field for the Polyakov loop. 
Such a consideration implies restoration of deconfinement when there is unbroken chiral 
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symmetry. This would lead to the Polyakov loop taking large expectation values inside 
the core of the string. Baryons are very heavy in the confined phase and very light in the 
deconfined phase. So baryons in the medium passed by a moving string would like to remain 
inside the string. This will give rise to increasing baryon number in the whole string loop as 
it collapses and moves through the medium. In the extreme case, the baryon number inside 
the string may become large enough to prevent the collapse of the string loop. The only 
way it can then decay is by quantum fluctuations involving the domain wall. 

There can be interesting effects if a defect network forms in the chiral transition in heavy- 
ion collisions. Since in the core of the string chiral symmetry as well as deconfinement are 
restored, it will have large energy density of the order of ~ 500 MeV/ fm 3 . For non-zero 
quark masses all the phase variation around the string will be confined to a domain wall 
bounded by the string. Similar domain walls in QCD with large chemical potentials have 



been discussed before |13j. Across the domain wall the U(1)a phase varies by 2tt. The 
surface tension of such a wall can be ~ 500Mev/ fm 2 , considering that its size is governed 
by the rj mass. So a string loop of radius 3/m will have a total energy ~ lOGeV. In 
case the dynamics of shrinking of the string loop happens after the freeze out, then the 
shrinking will be less dissipative. Finally the string loops will shrink to a region of lfm size 
depositing ~ lOGeV energy making it a hot spot. However an idea of size distribution of 
the string loops and their number can only be checked by real time lattice simulations of the 
linear sigma model. Such hot spots can give rise to a non-isotropic/non-azimuthal energy 
distribution. It will be interesting to see whether the non-isotropy due to this effect can be 
larger than the usual statistical fluctuations. 

In summary, we have argued that the symmetry breaking U(Nf) x U(Nf) — > U(Nf) 
allows for the existence of topological (abelian and non-abelian) string defects. We have 
studied the static configuration of these defects in the chiral sigma model. These defects 
may form during the chiral transition in the high energy heavy-ion collisions as well as in 
the early universe. We speculate that formation and subsequent evolution of the network 
of these string defects can give rise to inhomogeneous distribution of baryons and also the 
energy density. 
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